Peeling and Sliding in Nucleosome Repositioning 
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We investigate the mechanisms of histone shding and detachment with a stochastic model that 
couples thermally-induced, passive histone sliding with active motor-driven histone unwrapping. 
Analysis of a passive loop or twist defect-mediated histone sliding mechanism shows that diffusional 
sliding is enhanced as larger portions of the DNA is peeled off the histone. The mean times to histone 
detachment and the mean distance traveled by the motor complex prior to histone detachment 
are computed as functions of the intrinsic speed of the motor. Fast motors preferentially induce 
detachment over sliding. However, for a fixed motor speed, increasing the histone-DNA affinity (and 
thereby decreasing the passive sliding rate) increases the mean distance traveled by the motor. 

PACS numbers: 87.16. Sr,87.10.-l-e,87.16.Nn,87.14.Gg 



The nucleosome is comprised of double-stranded DNA 
that wraps l| turns around the edge of disk-shaped 
histone proteins. These structures impart a DNA 
accessibility-based code in addition to the sequence-based 
genetic code |T] . Positioning of the histone particles de- 
pends on the local histone-DNA affinity (via e. g., DNA 
bendability [1]) and can modulate the accessibility of 
DNA. The remodeling of nucleosomes, involving perhaps 
histone sliding and/or detachment occurs during cellular 
processes that require DNA accessibility [2j. For exam- 
ple, in replication, translation of a replication fork along 
the DNA requires histones be either shifted or removed 
[21 S] . It has been proposed that these processes may be 
passive, in which parts of the histone-DNA complex ther- 
mally unwind and rewind, leaving transient regions of the 
DNA accessible El El [7] . Histone sliding has also 
been observed and modeled theoretically [8j. However, 
since the histone-dsDNA binding energy is ~ iOkBT, 
sliding mechanisms must exploit low energy thermal ex- 
citation such as the propagation of extra-length loops 
or twist defects [S]. Although these thermal mecha- 
nisms result in slow histone sliding, more rapid nucle- 
osome remodeling can be catalyzed by molecular motor 
enzymes [TDl E] • A number of ATP-dependent remod- 
eling factors, such as SWI/SNF, ISWI, and SWI2/SNF2 
have been identified to be important during transcription 
[H El [H]. The action of these remodeling factors may 
be required for more rapid detachment of histones from 
the DNA substrate. Apart from the observation that 
chromatin remodeling factors have conserved helicase- 
like ATPase domains [T2], and that histones can both 
slide along and dissociate from DNA, there has been no 
explicit mechanistic hypothesis for enzymatically medi- 
ated chromatin remodeling. 

In this Letter, we propose and explore the conse- 
quences of a simple, explicit mechanism for enzyme- 
mediated chromatin remodeling. The remodeling factor 
is modeled as a processing motor that runs into an iso- 
lated dsDNA-wrapped histone particle. Such a motor 
not only can bias histone sliding, inducing a drift, it can 
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FIG. 1: The kinetics of single, isolated histone motion me- 
diated by thermal creation and propagation of loops. Flaps 
unbind and bind at either end of the contact region with rate 
fc„ and kfj, respectively. Closure of a flap into a loop occurs at 
rate a and loop ejection into a flap occurs with rate /?. Unbi- 
ased loop hopping rates in the interior of the contact region 
occurs with rate k. The asterisk labels a fixed site on the 
DNA substrate. A motor complex is shown to the left. 



also wedge itself underneath the histone, peeling off the 
DNA. The full stochastic process is described by an effec- 
tive master equation, from which we find the mean times 
to histone detachment, and the mean travel distance of 
the motor before it peels the histone off the DNA. Figure 
[Tj depicts the histone-DNA contact region, or footprint, 
and a thermal inchworming mechanism. Small flaps (or 
defects) are thermally excited with rate fc„ by fluctuat- 
ing segments of DNA that momentarily unwrap from the 
either edge of the histone-DNA contact region. The flap 
can rebind to the segment from which it detached with 
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rate or, it can reattach to the next-nearest neighbor 
segment with rate a, generating a loop. Once an interior 
loop is formed, it can hop to the right or left with rate 
k. The histone has translated one loop arclength only 
when a loop entering on one side of the contact region 
has left the opposite end. In this reptation-like mecha- 
nism, the net motion of an isolated histone is related to 
the statistics of diffusing loop "particles" that exit the 
end opposite from which they were thermally excited, 
analogous to the gambler's ruin problem [T3]. The sites 
of the coarse-grained lattice represent individual histone- 
DNA contact "bond" segments, separated by a distance 
equivalent to the typical arclength of a flap, about 10 
base pairs 

For a histone particle wrapped with n segments con- 
taining a single interior loop, we denote the probability 
that the loop is at position i {0 < i < n) hy Vi- Since the 
injection of loops is slow (a <C fcf,) [HlIIl], the probability 
for multiple loops is low, and loop-loop interactions can 
be neglected. The master equation governing the interior 
diffusive motion of the loop is Vi = kVi^i + kVi+i — 2kVi 
for 2 < i < n — 1. Near the left edge the bound- 
ary conditions are Vi = fcp2 — {P + k)V\ + aV^ and 
— kuVo — (fcfc + 0)7^0 -I- pVi . Similar expressions hold 
at the right edge. The probability that a loop exits the 
opposite end from which it entered can be found by con- 
sidering a unit source of probability Vq = 1, and solving 
the resulting recursion equations |13j . The total proba- 
bility of starting in state Vq and arriving in state Vn is 
{n-i)i3-^C+l3k(c.+2k,) - Since flaps are thermally generated 
with rate and produce the state P* with rate PVn, 
the effective hopping rate of the entire histone along the 
DNA is 



Pn 



akku 



{n - l)pkb + k{a + 2kb) 



(1) 



For an isolated, fully-wrapped histone, flap generation 
at each end is statistically independent and its motion 
is that of unbiased diffusion with diffusivity p^r. The 
loop hopping rate k will be approximated as two sequen- 
tial flap detachment and reattachment steps, without any 
additional barriers, giving k ~ k^kb / [k^ -\- k^) . Similarly, 
(3 can be approximated by the detachment rate. Finally, 
the insertion of a loop is proportional to the flap reat- 
tachment rate, but reduced by a factor associated with 
the bending required to form a loop. Thus, we use the 
physical rate estimates 



awQ!oA:fc, P ~ ku, and fc w fc„/ct/(fc„ -I- fcb), (2) 

where ao ~ e^^b/ksT ^ ^ scales the loop injection rate 
by the energy Eb required to bend the DNA into a small 
loop. These approximations make our results depend on 
just w/a and K = k^/ku, and the maximum number N 



of DNA-histone contacts. Since Eb 23kBT E] and 
ao ~ 10^^*^, we present all results in terms of quantities 
rescaled to the slowest rate a in the problem. 

Figure [2^ shows the normalized effective hopping rate 
Pn/oL as a function oi K = kb/k^ and n. For flxed kb, 
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FIG. 2: (a) Effective diffusion coefficient Pn/a as a function of 
flap binding constant K — kt/ku- (b) Diffusivity as a function 
of the degree of histone-DNA wrapping n. The parameters 
used are a — 10~^''fcb,/3 = and k — kbku/{kb + ku) (Eq. 
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Pn decreases with decreasing flap unbinding rate The 
effective hopping rate also depends on the number n of 
histone-DNA contact elements, as shown in Fig. |2]d. The 
higher the degree of wrapping of a histone, the larger the 
71, and the smaller its loop-mediated diffusion coefficient. 

Now consider a processing motor complex that moves 
unidirectionally with rate w along the DNA. For simplic- 
ity, we assume that the motor moves by steps of length 
approximately equal to that of a thermally generated 
loop. Fig. [1] shows a motor to the left of a histone- 
DNA contact region. The motor can advance only if the 
segment in front of it is cleared and not attached to the 
histone. The edge of the contact region just ahead of a 
motor can be cleared by a thermally-excited flap, allow- 
ing the motor to possibly move forward and slip under 
it. Consequently, the effective forward hopping rate w of 
a motor that is peeling a histone is reduced from the free 
hopping rate w by a factor dependent on the reattach- 
ment rate fcf,: 



w 



wku 
kh 



w - 



(3) 



The motor can also advance into the contact region by 
loops generated at the opposite end, but that have arrived 
after diffusing across the footprint. Thermal excitation 
of a flap at the far end, and subsequent propagation of 
a loop to the opposite end, occurring with rate p„, mo- 
mentarily creates a flap of two bond segments near the 
motor. Again, flap closing competes with motor inser- 
tion underneath the flap. As shown in Fig. |3] there are 
three fates available to a two-bond flap nearest the mo- 
tor. The motor can remain in its position and both bonds 
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FIG. 3: Details of the motor-flap interaction. When a loop 
generated on the opposite end exits the side nearest the motor 
(with efi'ective rate p„), the motor can move forward 0,1, or 
2 steps with relative probabilities /o,/i, and /2, respectively. 



of the flap can close (with effective rate kb/2), the mo- 
tor can move one position forward and one bond of the 
flap can close (with effective rate wkt/{w + kb)), or, the 
motor can move forward two sites (with effective rate 
w/2) under the flap, preventing it from closing at all. 
The total transition rate available to the two-bond flap 
is rx = kh/2 + kf,w/{kf, + w) + w/2, from which we find 
the probabilities for each pathway 



2rT' 



kbW 



and /2 = 



2rr, 



(4) 



Note that /o + /i + /2 = 1 and that the ratios of these 
pathways are determined by the motor speed w. These 
probabilities apply only to flaps on the motor side of 
the contact region. On the opposite, motor-free end, we 
assume that single bond flaps close sufficiently fast {K = 
kb/ku ^ 1) such that spontaneous desorption of histones 
is negligible in the absence of a processing motor. 

As the progressing motor partially unwraps the his- 
tone, the length n of the contact region decreases, in- 
creasing the effective sliding rate p„ of the remaining 
contact region (cf. Fig. |2|d). The sliding mechanism 
becomes increasingly important as the degree of unpeel- 
ing increases, and the histone becomes more efficient at 
diffusively escaping the motor. If m is the position of 
the motor, and n is its distance to the far end of the 
footprint, a master equation that incorporates all of the 
kinetics described above can be written for the probabil- 
ity P{m, n, t): 



P{m, n, t) = pm [P(m, n + 1) - 2P(m, n) + P(m, n - 1)] + w [P{m - 1, n + 1) - P{m, n)] , n> N + 1 

P{m,N,t) = -{w + pN)Pim, N) + wP{m - 1, N + 1) + foPN~lP{m, N - 1) + /ipivP(m - 1, N) + pN P{m, N + 1), 

P(m,n,t) = —(€) + pn)P(m,n) + wP{m - l,n + 1) + foPn-iP{m, n - 1) + fiPnPim — 1, n) 

+/2Pn+iP(m-2,n + l), 3<n<Af, 



(5) 



r 



with F(m, 1) = -kuP{m, 1) + wF(rn- 1, 2) + /2P2-P(to- 
2,2) and P(m, 2) = -{w + p2)P{m,2) + wP{m ~ 1,3) + 
fiP2P{,m-l, 2)+f2P3P{m-2, 3). The statistics of the hi- 
stone detachment time can be found by considering only 
the relative distance n, and when it flrst reaches n — 0. 
Upon summing the time-integrated probabilities over all 
possible motor positions m > 0, Eqs. [5] are succinctly 
expressed as 



fln+lQn+l — {o-n + bn)Qn + &n-lQn-l — —^n.N (6) 

where Q„ ee Jg^ dt X]m=o ^("^' ^^"^ = 
Sm=o 0) is the initial condition with the motor 

at the left end of the contact region of a fully wrapped 
histone. In Eq. [6] the transition rates are 



a„ = w + Pat, pn n> N 



an^W + f2Pn, bn = foPn 2 <n< N. 



(7) 



Conditions just prior to detachment {n ~ 1) also require 
bo = bi = and ai = fc„. 

The mean flrst detachment time of the process de- 
scribed by Eq. |6]is found by solving the corresponding 
Backward Kolmogorov equation [13], or by directly time- 
integrating the solution {T^} — X^^i Io° Qn(i)dt: 



'>n+k 



n=l 1=1 k=0 



(8) 



The mean distance (m^) traveled by the motor before 
it detaches the histone can be constructed from the time- 
integrated quantity i?„ = di ^^^q TOP(r7i, ti, t), 
which obeys an equation similar to (|6|: a„+ii?„+i — 
(o„ -I- bn)Rn + bn-iRn~i — Gn- The additional inho- 
mogeneity terms G„>jv = —Qn, Gn — ~wQn+i, and 
Gn<N = -{w + 2/2P«+i)Q„+i - fiPnQn, rely on the 
solution Qn to Eq. [6] Upon solving for we flnd 
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poo 

(md-l) = fc„/ dt'^mP{m,l,t) = kuRi. (9) 

In Fig. wc plot the mean first detachment times as 
a function of the motor speed, for A'^ = 15 and various 
flap binding constants K = kij/k^- For slow motors, the 
double flaps near the motor nearly always close (/q — s- 1) 
before the motor can advance. The histone unit diffuses 
along the DNA with rate pm and is only occasionally 
blocked by a relatively stationary motor. It spends most 
of its time away from the motor, preventing it from slip- 
ping under flaps that momentarily open. The resulting 
mean first detachment time diverges exponentially in the 
w/a — > limit. 
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FIG. 4: (a) Tlie normalized mean time to detacliment, 
a{Td)/N, as a function of tlie normalized motor speed w/a, 
for various K. The dotted curve shows the dramatic decrease 
in detachment times corresponding to isT = 20 and a hypo- 
thetically fixed small difFusivity pjv that suppresses sliding, 
(b) The mean distance {rrid} traveled by the motor before it 
peels off the histone. In both plots, the maximum wrapping 
length was set to N = 15, corresponding to a flap/loop length 
of about 10 base pairs. Since we have assumed tight binding 
(large K), and ignored spontaneous histone desorption, our 
results (Td) and (mj) for smaller values of K are only upper 
bounds. 

In the large w/a limit, the remodeling complex slips 
under the flaps before loops can propagate through the 
contact region and the peeling mechanism dominates. As 
w is increased, the peeling time decreases. In the ex- 
tremely fast limit w/kh 3> 1, the motor nearly always 
inserts itself underneath flaps before they can reattach 
(/o ~ /i ~ 0, and /2 ~ 1) and Eq. |8] simplifies to 

(Td) = + En=2(fc^' + Pn)~^ < Nk-\ The mean 
detachment time of a histone subjected to an infinitely 
fast motor would be (T^) = N/ k^ if flaps originated only 
in front of the motor. However, flap generation at, and 
loop propagation from the other end of the contact re- 
gion enhances the arrival of flaps near the motor, thereby 
decreasing {Td) to values below N/ku- 

In Fig. |4Jd, we plot the mean distance traveled by the 
motor before it detaches a histone, as a function of w/a 
for various K . For low motor speeds, sliding dominates 



peeling, and the motor moves significantly. Note as the 
fiap binding affinity K increases, the mean distance trav- 
eled by the motor also increases. For flxed k^, increasing 
K by decreasing fc„ leads to a dramatic increase in (T^) , 
especially at low motor speeds. This extra time allows 
the histone to slide a longer distance before being peeled 
off by the motor. Although increasing K reduces the 
effective sliding rate of an isolated histone, mechanisti- 
cally, the increase in mean distance traveled arises from 
the relative suppression of double flaps near the motor, 
slowing down the peeling process. 

For larger w/a, peeling becomes dominant and the mo- 
tor moves a mean distance {nid) ~ N before scrapping 
the entire histone off. Again, in the w/fcb 2> 1 limit, 
an analytic form can be explicitly found: {md) /N — 

l + X]^=2Pn/(^"~'"f^")- This slow increase in the normal- 
ized mean distance traveled, {rrid) /N, with N arises from 
the relative reduction of double flaps near the motor, 
thereby diminishing the frequency of this efficient mode 
of peeling. This reduction occurs for larger N through 
the loop propagation probability p„. 

In summary, we have presented and analyzed a 
stochastic model for ATP-dependent chromatin remodel- 
ing that incorporates a competition between sliding and 
peeling of a single isolated histone. Our main findings in- 
clude formulae for the mean detachment time and travel 
distance, as functions of motor speed and binding affin- 
ity. Histones with larger binding affinity to the DNA 
substrate slide farther before being peeled off by the mo- 
tor. We have neglected effects such as DNA sequence 
dependence (which might affect w,k^,kh, and k) and 
twist defect mobility [5]. This latter mechanism may 
only open much smaller end flaps that allow motor in- 
sertion and peeling. Nonetheless, our results are useful 
for dissecting and quantifying in vitro measurements of 
passive sliding and ATPase-assisted nucleosome remod- 
eling, particularly when physical chemical conditions can 
be controlled and w, kj, tuned. 

Nucleosomes in vivo typically consists of an array of 
interacting histones [TTJ[TS]. Thus, our results would be 
valid only if the typical DNA linker length £ between 
two adjacent histones were larger than {m,d) + N . If 
(- < {rud) + N, the histones would be squeezed together 
by the motor. In such cases, peeling would be enhanced 
as the sliding mobility is reduced due to histone-histone 
exclusion. Finally, it should be mentioned that the anal- 
ysis of Eq.[5] is directly applicable to the study of dehy- 
bridization and driven translocation of double-stranded, 
mononucleotide nucleic acids through nanopores |16j . 
Our results, after reidentification of the parameters, pro- 
vide an analytic solution to the mean passage times of 
translocation in the high driving force limit. 

This work was supported by the NSF (DMS-0349195) 
and by the NIH (K25AI41935). 
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